Covariance matrices and the separability problem 
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We propose a unifying approach to the separability problem using covariance matrices of locally 
measurable observables. From a practical point of view, our approach leads to strong entanglement 
criteria that allow to detect the entanglement of many bound entangled states in higher dimensions 
and which are at the same time necessary and sufficient for two qubits. From a fundamental perspec- 
tive, our approach leads to insights into the relations between several known entanglement criteria 
- such as the computable cross norm and local uncertainty criteria - as well as their limitations. 
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Entanglement plays a central role in applications of 
quantum information science as well as in the foundations 
of quantum theory. Quite naturally, one of the problems 
that have received a significant amount of attention is 
the question to decide whether a given state is entangled 
or separable. In fact, the development of separability cri- 
teria 0, 0, H, 0, @, Q has been one of the key activities in 
quantum information theory: On top of certifying a given 
state to be entangled, often in an experimental context, 
they often provide physically valuable intuition concern- 
ing the structure of the entanglement in a given state. 
Among these attempts^methods using positive maps [l[ 
or convex geometry [1, H, 0] turned out to be fruitful. 

In this work, we propose a unifying approach to find- 
ing criteria for separability for finite dimensional systems 
in terms of covariance matrices (CMs). In the infinite- 
dimensional setting (in particular for Gaussian states) 
such CMs constitute a well-established and powerful tool, 
not least due to the experimental accessibility of quadra- 
ture measurements using homodyning @, B Q • In con- 
trast, for finite-dimensional systems, the theory is yet 
hardly developed [nj 11 \ ■ We introduce a framework for 



CMs for finite-dimensional systems, formulate a general 
separability criterion and evaluate it for various scenar- 
ios. The merits of this approach are two-fold: (i) The re- 
sulting criteria are very strong and allow, notably, to de- 
tect many bound entangled states, (ii) Our approach pro- 
vides a framework to link and understand several existing 
criteria like a recent criterion using the Bloch represen- 
tation 11211 . the one based on local uncertainty relations 
(LURs) [13], or the cross-norm or realignment (CCNR) 
criterion [2|, the latter being an immediate corollary of 
our theory. 

The main idea. - Let us start by defining CMs. Let g 
be a given quantum state and let {Mk : k = 1, . . . , N} be 
some observables. Then the N x N CM 7 - dependent 



on the state g and the choice for {Mk} - is given by 

7id = ((MiMj) + (M 3 -Mi»/2 - (Mi)(Mj). (1) 

This is a real, positive definite matrix [l4j and its diago- 
nal entries are just the familiar variances, 7^ = 5 2 {Mi) e . 
The CM has a useful concavity property: if g = J2kPkQk 
is a convex combination of arbitrary states, then 



7(0) > ^Pkiigk)- 



(2) 



This can be shown to emerge from the generating func- 
tion of the moments @, [H| > reflecting the fact that the 
variance of an observable increases under mixing. 

Now, let g be a bipartite state on Tt = Ha®Hb, where 
d,A {d B ) is the dimension of Ha (Hb). We can choose d\ 
observables {Ak} on Ha such that they form a Hilbert- 
Schmidt orthonormal basis of the observable space, that 
is, they obey TrL4fcv4;] = Sk.i- E.g. for a qubit, we may 
choose normalized Pauli matrices including the identity. 
Similarly, we may take {B k } as a basis of observables in 
Hb, and consider the total set {M k } = {A k ®l, l(g)Bk}. 
The CM 7 then has the block structure 



lie, {Mk}) 



A C 
C T B 



(3) 



where A = j{g A ,{A k }) and B = j(g B ,{B k }) are CMs 
of the reduced states, and C has the entries Ci.j = (Ai <g> 
Bj) — (Ai)(Bj). This matrix will form the starting point 
to characterize the separability properties of g. 

To formulate the separability criterion, recall that a 
state is separable iff is a convex combination of prod- 
uct states, i.e., g = J2kPk\ a k,bk)(ak,b k \. We then use 
Eq. ([2]) and the fact that for a product state the block C 
in Eq. ((3|) vanishes, to arrive at the following observation: 
Observation 1 (CM criterion). Let •y(g) be a CM 
as in Eq. f3|). If Q * s separable, then there exist states 
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\a k )(a k \ on Ha, \h)(bk\ 
Pk such that for ka = 
Y.kPkl(\bk)(bk\) we have 



on H.b, and convex weights 
Pkl{\ak){a k \) and k b = 

'1% 



7(e. Wk}) >k a ® 



kb- 



(4) 



If no such ka/b exist, g must be entangled. 

We refer to this criterion as covariance matrix criterion 
(CMC) [lij]. Subsequently, we will show that this con- 
dition can be made an efficient and physically plausible 
test. We will see that as such the criterion is independent 
of the choice of the observables {Ak}, {Bk}, but that a 
certification of a violation of the criteria is simplified by 
choosing the Schmidt-basis in operator space or by first 
using an appropriate local filtering. 

Let us first note some facts concerning the CM 7 
and the detailed presentation will be given 

elsewhere [l7| . A general change of the observables 
{M k } i-> {M k } with M k = J2iVk,iMi gives rise to a 
map 7(0, {M k }) = ^7(0, {M k })fi T . From this it is easy 
to see that the CMC is, in principle, independent of the 
choice of the orthonormal {A k }, {B k } , however, a suit- 
able choice of them simplifies the falsification of Eq. ((3]). 
A unitary transformation g 1— ► UgU* induces a trans- 
formation M k 1 ► Mfc = WM k U = J2i°kdMi where O 
is orthogonal. The orthogonality of O is equivalent to 
the statement that the {M k } are also orthogonal So 
the eigenvalues of 7(0, {M k }) are invariant under unitary 
transformations of the state. This also means, notably, 
that by a suitable choice of the observables, the matrix 
C in Eq. Q can be made diagonal by a singular value 
decomposition. If g is a rf-dimensional pure state and the 
{Mk] are orthogonal, then 7 = P/2 where P = P 2 is a 
projector onto a 2{d — l)-dim subspace of the total d 2 - 
dim space. This can be directly calculated for a specific 
state and some {M k } (HI- Then, the general statement 
follows from the second property. This also implies for 
Eq. (J) that Tt[ka] = d A - 1 and Tt[k b ] = d B - 1. 

Evaluation for two qubits. - As a first example, we 
take as observables the normalized Pauli matrices, i.e., 
{Ak} = {B k } = {l,a x ,a y ,a z }/V^. Then, in the def- 
inition of 7 in Eq. ([3]) two rows and columns (corre- 
sponding to the A\,B\ = l/v2) equal zero, thus it suf- 
fices to consider 7 as a 6 x 6 matrix, originating only 
from {Ak} and {Bk} for fc = 2,3,4. To characterize the 
Ka, note that for a pure state on Ha the 3x3 ma- 
trix "f(\a) (a\, {Ak}) is a two-dimensional projector, i.e., 
7 = P/2 = (1 3 - \<f> a )(<f> a \)/2, where |0 O ) e R 3 . In turn, 
each matrix of the form (I3 — \<p)(<p\)/2 is a valid CM, 
since any such matrix can be obtained from a special 
(I3 — |</> a )(</> Q |)/2 by an orthogonal transformation (in 
R 3 ), and for the special case of a qubit there is a one to 
one correspondence between such orthogonal transforma- 
tions and unitary transformations of the state [3] . From 
this, it is easy to see that K A = (I3 — p A )/2 where pa is 
a real density matrix on R 3 , and we can summarize: 



Proposition 2 (Qubit criterion). Let g be a state of 
two qubits, and 7 be the 6x6 CM as in Eq. {3J) with 
{Ak} = {Bk} = {a x ,a y ,a z }/V2- Q fulfills the CMC iff 
there exist 3x3 density matrices pA, Pb with 



1 -1q/2+{ Pa ®pb)/2>Q. 



(5) 



If we find complex pa/b-, their real part saturates Eq. (J5]) 
as well. Finding the pa/b is a simply solvable semi- 
definite feasibility problem Eq. (O will become im- 
portant later in the context of LURs. But let us discuss 
the general case first. 

Evaluation of the CMC for the general case. - Let us 
first assume that dA = ds — d. 

Proposition 3 (General criterion). Let g be a state 
with dA — ds = d and A, B, C be as in Eq. |3|). If g is 
separable, then 



,1- 



2> '\Ci,i\ < (Tr[A] 



+ l) + (Tr[S] -d+1) 
(l-Tr[^]) + (l-Tr[e 2 B ]). (6) 



Proof. First note that a necessary condition for a 2 x 2 
a c 



matrix X 



to be positive is that 2\c\ < a + b. If g 



were separable, then by the CMC we have Y = 7 — ka ® 
kb > 0. This implies that all 2 x 2 minor submatrices 
of Y have to be positive. Hence, for all i,j 2\Cij\ < 
Ai i + Bj j — (nA)i,i ~ ( K B)j.j ■ Summing over all i, j proves 
the claim, using Tr[« A ] = Tt[kb] = d - 1. Eq. © 
always holds since X^fc^fe = dl [21(, hence A4 j = 
Ekt(Ak) 2 =d-Tr[g 2 A ]. □ 
As explained before, one may choose C diagonal, so 
considering the diagonal only is no restriction of gener- 
ality (2ol |. One may further try to improve Proposition 
3 by taking 4x4 submatrices into account [T3|- Physi- 
cally, Proposition 3 states that if the correlations Ca are 
large, then g is entangled. The question remains how to 
find those observables for which the C^, are large. This 
problem can be overcome by making use of the Schmidt 
decomposition in operator space [2l| . Recall that we can 
express any state as g = J2k ^fc(^fe ® ^k )• Here, Afe > 
for all fc, and the {Gk} form an orthonormal basis of the 
operator space. Denoting g^ B = Tr[G^ S ] we can now 

use the set {G^ B } in Proposition 3: 

Proposition 4 (CMC for Schmidt form states). If 

g is separable, then 



2£|A fe 

fc 



A 2 fc5fcVl 



< 2 



A 



This is a direct application of Proposition 3 to the diag- 
onal of C. One can find examples of states, which are de- 
tected by the CMC via Proposition 2, but not by Propo- 
sition 4 [22|]. However, Proposition 4 connects now the 
CMC to the CCNR criterion @,[2l|: 
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Corollary 5 (Connection to the CCNR criterion). 

If g is separable, then Afc < 1. 

This follows from Proposition 4 and the general re- 
lation a 2 + b 2 > lab. Since the condition J^k < 1 is 
just the CCNR criterion, this implies that the well-known 
CCNR criterion is a direct corollary of our theory [23j |. 

Enhancing the criterion via local filtering. - A general 
strategy strengthening the presented test is the joint use 
of the CMC with local filtering operations. Such opera- 
tions are maps of the form g i— > g = (F a ®Fb)q(F a ®Fb)'' ! 
where the F A / B are arbitrary invertible matrices. They 
preserve the entanglement or separability of a given state. 
Local filtering transforms a generic (full rank) g into 



FNF is separable, then 



d A d B 



d\-l 

£ 

=i 



1 + Y] UGf ® Gf) 



(7) 



where & > for all i and the {G^ B } are traceless 
and orthogonal observables 0, O, I2II. The matrices 
Fa,Fb can be found constructively [27[. We will re- 
fer to this form as the filter normal form (FNF). The 
extraordinarily helpful property is that for a state in 
the FNF with the {Gf^ B } as observables the blocks 
in the CM are diagonal. We have then in Eq. ([3]) 
A = diag(0, 1,1,..., l)/d A , B = diag(0, 1,1,..., l)/d B , 
and C = diag(0, £1, £2, ■ • ■ , £d 2 -lVO^^s), and obtain: 

Proposition 6 (CMC under filtering). If a generic 
state g is separable and dA = d B = d we have in its FNF 
Yt\ l Zi<d*-d. 

This is a very strong criterion for separability, as some 
examples will show below. Interestingly, it is also nec- 
essary and sufficient for two qubits [281 ] : For them, it is 
easy to see that the {Gf^ B } in Eq. ([7]) are effectively the 
Pauli matrices, and the separable states of this form are 
known to be an octahedron inside some tetrahedron [l9| , 
the borders of which are described by Proposition 6. Al- 
ternatively, this can be seen from the results in Ref. [2(| . 

Let us consider the asymmetric case, when d A < d B - 
The CM 7 in the FNF is then similar as before, however A 
and B arc not of the same dimension. Two observations 
arc now helpful: First, we do not sum over all On the 
other hand, we cannot subtract all of the (K B )i,i anymore, 
since d 2 B — d\ diagonal elements of k b do not occur in the 
sum. Moreover, 7 consists of a linear part 7.^ = (M{Mj + 
MjMi)/2 and a nonlinear part 7^ = {M i ){M j ). Since 
the linear part is compatible with convex combinations, 
the linear part on the right hand side of Eqs. ([21 14|) has 
to coincide with the linear part on the left hand side. 
The non-vanishing elements of 7 origin from the linear 
part only, so the missing {n B )i^ are {n B )i,i = l/d B — 
J2kPk(Gi) 2 bk ^ bk \, hence (K B )j,j < l/d B . So, if g in the 



6 < d A d B {l-l/d A + (d 2 A-l)/d B (8) 

+ min{0; -(da - 1) + (d 2 B - d\)/d B })/2. 

The possibility of taking zero in the minimization comes 
from the fact that one may also omit the summation over 
[n B )is. Especially if d A -C d B this yields better bounds. 

It is instructive to compare this with the CCNR cri- 
terion which requires for separable states in the FNF 
< d A d B — {d A d B y/ 2 , and the dV-criterion [13], 

1/2 

requiring < (d A d B (d A - l)(d B - 1)) ■ For d A = 

d B all three criteria coincide. One can directly see that 
for dA < d B Eq. (|8|) is stronger than the CCNR crite- 
rion, but the dV-criterion is also stronger than CCNR. If 
d B — dA is small, the dV-criterion is slightly better than 
Eq. JSJ), for dA <C d B , however, Eq. ([8]) is drastically 
better than the CCNR and dV-criterion. 

Connection with the L URs. - LURs allow detection of 
entanglement in the following way [l3| : One takes observ- 
ables {Ak} and {Bk} on Alice's (resp. Bob's) space and 
computes two bounds Ua,U b such that X^fc^ 2 (^<0 — 
Ua, S/c^ 2 (-^ fc ) — Ub- Then for separable states 
J2 k S 2 (A k <g> 1 + 1 <g> B k ) > U A + U B holds. Physically, 
LURs show that separable states inherit the uncertainty 
relations from their reduced states, which is not the case 
for entangled states. For a given state g, however, it is 
usually not clear which {Ak} and {Bk} are suitable to 
detect its entanglement. We can formulate: 
Proposition 7 (Connection to LURs). A state g 
violates the CMC iff it can be detected by a LUR. 

Proof. The proof is given in the Appendix. □ 
This equivalence has two consequences: Results con- 
cerning LURs can be transferred to the CMC. E.g., 



Ref. [2l| gives now an alternative proof that the CMC 
is stronger than the CCNR criterion. Also, our results 
concerning the CMC allow to gain new insights in LURs: 
Corollary 8 (LURs for two qubits). There exist en- 
tangled two-qubit states which can not be detected by a 
LUR. 

Examples of such states can be found numerically [13] , 
showing that the filtering indeed improves the CMC. 

Examples and extensions. - Firstly, we take the 3x3 
bound entangled states arising from an unextendible 
product basis, mixed with white noise [29| . These states 
Qupb(p) are detected by Proposition 6 for p > 0.8723 
while the best known positive map detects them only for 
p > 0.8744 [21]]. Secondly, we checked randomly gener- 
ated 3x3 chessboard states [3fj| . The CCNR criterion 
detected 18.1 % of them, Proposition 4 detected 19.1 % 
and Proposition 6 detected 97.8%. 

Finally, let us note that the theory developed in 
this paper may be complemented by considering non- 
symmetric CMs [l7|- That is, one may define 7,^ — 
(MiMj) — (Mi)(Mj). For such matrices (which are now 



4 



hcrmitian, but still positive) one can directly derive a 
statement corresponding to Proposition 1. This criterion 
implies then the CMC as given for a symmetric 7. 

Conclusion. - We have proposed to investigate the sep- 
arability of finite dimensional quantum states using the 
CM for certain observables. We have demonstrated that 
this approach reveals the entanglement of many states 
and can lead to new insights into already existing crite- 
ria. It hence provides a further systematic framework of 
studying entanglement for composite quantum systems. 

We thank H.J. Bricgel, J.I. Cirac, R. Horodccki, M. van 
den Nest, G. Toth and an anonymous referee for discus- 
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QAP), the QIP-IRC, Microsoft Research, the DFG (SFB 
407), and the EURYI for support. 

Appendix. - Here, we prove Proposition 7. Note 
that the CM can be used to compute variances: if 
N = J^k^kMk is a linear combination of the M k then 
S 2 (N) = EijvaiMJiM = (i/| 7 (M fc )|i/>. If q vio- 
lates the LURs we can find Ak and Bk as above. We 



E 



and B k = J2iPi k)B i with A 



1 "i 



write Ak 

and Bk as in the definition of 7 in Eq. (O, leading to 
S 2 (A k 811 + 104) = © /? (fe) |7|a (fc) © /3 (fc) ). For 

the k a © k b we have k a © k b ■= J2kPkl(\ a k) («fc|) © 
l(\bk)(bk\), hence (a^ fc ) © /3^\ka © «s|aW 
£ , Pl (5 2 (A k ) 1 ai ) (oi 1 + 5 2 {Bk ) 1 a! > (o| 1 ) . If the CMC were ful- 
filled, summing over k would yield J^k ^ 2 (^fc © 1 + 1 © 

Bk) > ElPl(Y,k S2 (Ak)\a k ){a k \ + Efc^ 2 (^fc)|6 fc >(fc fe | 



> 



min |a)(a|(Efc <5 {A k )\ a){a \) + mm lb){b \(J2k S ( B k)\b)(b\) > 
XJ a + Ub, which contradicts the violation of the LURs. 

Conversely, let us define X as the set of all matrices 
which can be written as ka © kb + P with some ka 
and Kb as in Observation 1 and a positive P. In other 
words, the CMC states that for separable states 7 G X. 
Geometrically, X is a closed convex cone. According to a 
Corollary of the Hahn-Banach Theorem for each 7 X 
there must be a symmetric matrix W and a number C 
such that Tr[W 7 ] < C while Tr[Wx] > C for all z G 
X. Since Tr [WP] > for all P > 0, we have W > 0. 
Now we use the spectral decomposition and write W = 
EfcAfclVfcXV-fcl =E fe A fc |a«ffi/3W}(aWffl/3( fc )|. Defining 



and Sfc = V%c"Ez Pi^ Bi, we have 
E fc <5 2 04fc©l' 



-4fc - \AfcE; "j 

for g that Tr[VF 7 ] = £ fc <^(4fe © 1 + 1 © B k ) < C. 

Furthermore, we have that all ka © kb G X and even 
more, we have that all 7a©7b G X . Hence, for a product 
state q = g A © g B we have C < Tr[PF( 7j 4 © 7 b)] = 

(E fc ^ 2 ' 
mii v (Efc' 52 ^) 



fcU (Afe) e A +Efc^ 2 (^fc)e s ) wn i cn implies that C < 



leading to a violation of the LURs. 



-min eB (j: k 6 2 (Bk) e s)=:UA + 



□ 
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